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ABSTRACT 

We investigate the impact of the observed correlation between a galaxies shape and its sur- 
rounding density field on the measurement of third order weak lensing shear statistics. Using 
numerical simulations, we estimate the systematic error contribution to a measurement of the 
third order moment of the aperture mass statistic (GGG) from three-point intrinsic ellipticity 
correlations (III), and the three-point coupling between the weak lensing shear experienced by 
distant galaxies and the shape of foreground galaxies (GGI and GII). We find that thkd order 
weak lensing statistics are typically more strongly contaminated by these physical systemat- 
ics compared to second order shear measurements, contaminating the measured three-point 
signal for moderately deep surveys with a median redshift ^ 0.7 by ^ 15%. It has been 
shown that accurate photometric redshifts will be crucial to correct for this effect, once a 
model and the redshift dependence of the effect can be accurately constrained. To this end 
we provide redshift-dependent fitting functions to our results and propose a new tool for the 
observational study of intrinsic galaxy alignments. For a shallow survey with ^ 0.4 we 
find III to be an order of magnitude larger than the expected cosmological GGG shear signal. 
Compared to the two-point intrinsic ellipticity correlation which is similar in amplitude to the 
two-point shear signal at these survey depths, third order statistics therefore offer a promising 
new way to constrain models of intrinsic galaxy alignments. Early shallow data from the next 
generation of very wide weak lensing surveys will be optimal for this type of study. 
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1 INTRODUCTION 

Weak gravitational lensing represents a powerful tool to investi- 
gate the large-scale distribution of matter. The majority of lensing 
results to date have focused on using two-point shear statistics to 
constrain the matter density parameter fim and the matter power 
spectrum normalisation as (van Waerbeke et al. 2001; Hoekstra et 
al. 2002; Bacon et al. 2003; Jarvis et al. 2003; Hamana et al. 2003; 
Rhodes et al. 2004; Heymans et al. 2005; van Waerbeke et al. 2005; 
Massey et al. 2005; Semboloni et al. 2006; Hoekstra et al. 2006; 
Massey et al. 2007; Benjamin et al. 2007; Fu et al. 2008). As these 
two parameters are strongly degenerate however there is great in- 
terest in measuring higher order statistics as their combination with 
the two-point shear statistics can effectively break the degeneracy 
between Q^n and as (Bernardeau et al. 1999). To date, there have 
been few measurements of three-point shear statistics (Bernardeau 
et al. 2002; Jarvis et al. 2004; Pen et al. 2003). These results found 
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that the three-point shear statistics were significantly affected by 
systematics even when the two-point statistics showed a very low 
systematic level. The aim of this paper is to investigate, by using 
ACDM N-body simulations, whether the intrinsic alignment of the 
sources (see for example Heavens et al. 2000) and the correlation 
between the shear field and the intrinsic ellipticity of the sources 
(Hirata & Seljak 2004; Heymans et al. 2006; Mandelbaum et al. 
2006; Hirata et al. 2007) can explain the presence of some system- 
atics. 

In the weak lensing regime the observed ellipticity Co of a 
source galaxy is related to the original ellipticity es through: 

Bo ~ Es + 7, (1) 

where 7 = 71 + i72 is the complex shear, and e is the complex 
ellipticity defined as 

l-^exp(2i0) 

I+/32 R ■■ ^' 

where (j) is the angle between the semi-major axis and the a;-axis, 
P is the ratio between the semi-major and semi-minor axis and R 
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is the response of a galaxy to weak lensing shear field. Note that in 
the commonly used KSB method (Kaiser et al. 1995) the respon- 
sivity R is expressed by the polarizabiUty tensor which is com- 
puted for each object and relates the measured weighted ellipticity 
to the shear field. In this analysis where no weight is applied we are 
able to calculate the average responsivity for our sample finding 
R = 0.89 (see Rhodes et al. 2000, Bernstein & Jarvis 2002). 

The two- and three-point ellipticity correlation functions be- 
tween the observed galaxies a, b and c are given by : 



(e>S) = {eX'> + GI+(7V>, 
{eoelel} = (ete'se's) + GGI + GII + (i^lS"), 

where the GI, GII and GGI terms are: 

GI = {etf") + {eh"), 

GII = (7"e^e=) + {j'etet) + {^^etel), 
GGI = {-n'el) + (7^=6?) + (7^7" e') • 



(3) 
(4) 

(5) 
(6) 

(7) 



One assumes that galaxies are intrinsically randomly oriented so 
the measurement of the observed ellipticity is an unbiased estima- 
tor of the shear. In a similar way one assumes that the measurement 
of the correlation between ellipticity of pairs and triplets of galax- 
ies is an unbiased estimator of the second and third order moments 
of the weak lensing shear field. However, this assumption is not 
correct for galaxies that are physically close as an intrinsic align- 
ment between the galaxies can be induced by a tidal field which 
acts on cluster scales as has been observed by Brown et al. (2002) 
and Mandelbaum et al. (2006). These results imply that the first 
term of Equation (3), hereafter referred as to 11, and the first term 
of Equation (4), hereafter referred as to III, are non-zero and they 
systematically affect the value of the two- and three-point weak 
lensing shear statistics estimated through the observed ellipticity 
alignment. It has been shown (King & Schneider 2002, 2003, Hey- 
mans & Heavens 2003, King 2005) that this bias can be removed 
if close pairs of galaxies are downweighted in the lensing analy- 
sis, something which is possible only when the redshift for each 
individual galaxy can be reliably estimated. For moderately deep 
surveys such as CFHTLS Wide (Fu et al. 2008), the intrinsic align- 
ment does not significantly affect the estimation of the normalisa- 
tion of power spectrum an . However, once the redshift information 
is added it is possible to exploit all the information carried by the 
shear signal, by using three-dimensional measurements. The mea- 
surement of the two-point shear statistics in redshift bins, i.e. to- 
mography (Semboloni et al. 2006) and the 3D-lensing (Heavens et 
al. 2006, Kitching et al. 2007), are particularly promising to con- 
strain the equation of state of dark energy, but are also more suscep- 
tible to stronger contamination by intrinsic alignment (Heymans et 
al. 2004, King 2005). 

A more subtle effect is the correlation between the induced 
weak lensing shear and the foreground intrinsic ellipticity field, also 
called "shear-shape" correlation (Hirata & Seljak 2004). This effect 
can be explained as follows. Consider a background galaxy with 
redshift a" whose light is deflected by a foreground over-density at 
. Consider now a galaxy at redshift stretched by the tidal forces 
generated by the dark matter over-density. Then the term (7°es) 7^ 
0. Similarly, for the three-point shear statistics, one can consider 
two sources with redshift and whose light is deflected by the 
same halo hosting a galaxy at redshift a'^ < z", z^: in this case the 
term (7°7''es> / 0. Finally, one can imagine a background galaxy 
with redshift being sheared by an over-density hosting galaxies 
at redshifts z'' and z" with the result (7"eseJ) 0. 

The first measurement of the shear-shape systematic effect 



on the Sloan Digital Sky Survey (SDSS) data by Mandelbaum et 
al. (2006) has been refined by Hirata et al. (2007), which uses 
subsamples from the SDSS and 2SLAQ surveys to check the de- 
pendence of the shear-shape correlation on the color, morphology 
and luminosity of galaxies. Hirata et al. (2007) predict that the 
shear-intrinsic alignment coupling could cause an underestimate of 
the normalisation of the power spectrum erg of up to ten percent. 
Heymans et al. 2006 (hereafter H06) provide a simple toy model 
to investigate this effect using N-body simulations which we de- 
scribe in more detail in section 3. This simple model also predicts 
a shear-intrinsic aligimient of 10% of the expected weak lensing 
shear signal, and provides a good fit to the measurements of Hi- 
rata et al. (2007). In contrast to the intrinsic alignment systematic 
that affects the lensing measured in the same tomographic redshift 
slice, the shear-shape systematic affects galaxies in very different 
tomographic redshift bins. The shear-shape systematic strongly re- 
duces the precision with which one is able to constrain the dark en- 
ergy equation of state by tomography, as shown by Bridle & King 
(2007). 

Adopting the same strategy as H06, we analyse a set of ACDM 
N-body simulations, in order to compare the amplitude of the in- 
trinsic ellipticity and the shear-shape terms to the three-point weak 
lensing shear statistics. 

The paper is organised as follows. In the section 2 we define 
the quantities and we describe the method used in this work to mea- 
sure three-point statistics. In section 3 we describe the simulations 
used in this work. In section 4 we present measurements of the 
three-point intrinsic alignment, and in section 5 we show the evolu- 
tion of the three-point coupling between shear and intrinsic elliptic- 
ity fields as a function of the redshift distribution of the lenses and 
sources. We also provide a fitting formula for the three-point shear- 
ellipticity correlation from different galaxy models. We conclude 
in section 6. 



2 THREE POINT SHEAR STATISTIC MEASUREMENT 

Following the approach of Pen et al. (2003) and Jarvis et al. 
(2004) we define the complex 'natural components' of the three- 
point shear correlation functions (Schneider & Lombardi 2003) 
for a triangle of vertex Xi, X2 and X3 and separations vectors 

ri = X3 - X2 ra = Xi - X3, ra = X2 - Xi: 

ro(ri,r2,r3) = (7(ri)7(r2)7(r3)e[-''(^i+^^+^^"), (8) 

ri(ri,r2,r3) = {7*(ri)7(r2)7(r3)e[-''f-^i+'^^+^=*'l>, (9) 

r2(ri, r2, rs) = {7(ri)7*(r2)7(r3)e[-''f*i-*=+'*=''l), (10) 

r3(ri, r2, rs) = (7(ri)7(r2)7*(r3)e[-''(*i+^^-*=''l), (11) 

where 7 is the complex shear and we indicate with 7* the com- 
plex conjugate of 7. The choice of the directions (f)i along which 
one projects the shear are free as long as the value of F, does not 
depend on the triangles orientation. Assuming that the Universe is 
homogeneous and isotropic, the eight real components of the shear 
correlation function depend on the side-lengths of the triangle, ri, 
r2, and rs. Thus we prefer to use a derived statistic which is eas- 
ier to visualize and interpret, namely the third order moment of the 
aperture mass statistic. 

The complex aperture mass is defined as : 



g(-2i«^ (12) 



Mi, = Map +iMx = J (frQ^{r)-/ir) 
where Q-&{r) is a filter of characteristic size -d. Using the aperture 
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mass statistics allows one to uniquely separate the E-mode (Map) 
and B-mode (Af x ) component of the measured shear (Crittenden et 
al. 2002), providing a powerful test for systematics. Indeed, weak 
lensing shear fields are E-type fields; thus the presence of other 
sources of distortion can be revealed by measuring Mx statistics 
which are non-zero only for B-mode fields. 

In practice, the variance and the third order moment of the 
aperture mass can be measured as a function of the two- and three- 
point shear statistics, respectively (Crittenden et al. 2002; Pen et 
al. 2002). Estimating the third order moment of the aperture mass 
through three-point correlation functions is preferred however, as 
it yields a higher signal-to-noise ratio than the direct measurement 
of the moments of the integral of Equation (12) (Schneider et al. 
2002). In addition, in the case of real data, where masked regions 
are present, the estimate of the aperture mass statistic through the 
measurement of the correlation function is unbiased, whilst the di- 
rect measurement of the moments is biased. By choosing the fol- 
lowing filter. 



Q4r) = 



2^^ J 



exp 



(13) 



the four components of the third order moment 
(M'ip), {MfpMx),(MapMx) and (Mj) can be obtained through 
integration of the three-point shear statistics Vi (see for example 
equations (45), (50) of Jarvis et al. 2004 and equations (61-71) of 
Schneider et al. 2005). There are other advantages of using aperture 
mass statistics defined by using Equation (13). This filter has infi- 
nite support but the exponential cutoff allows one to assume a finite 
support for real calculations. Furthermore, the Fourier transformed 
filter. 



2ti 4 



exp(-?7^) 



(14) 



were we defined 77 = si?, is a very narrow window filter, probing 
essentially modes with s ~ 2 /■&. 

We measure the four complex Ti components using a binary 
tree code built following the model suggested by Pen et al. (2003). 
We assign to each box an ellipticity (ei, 62) and a position (x,y) 
which are the average of the ellipticity and positions of the galaxies 
contained in the box. To each box we assign also a weight which is 
the sum of the weights assigned to each galaxy. The characteristic 
size li of each box is chosen to be the distance between the centre 
of the box and the furthest galaxy contained in the box. The cor- 
relation between triplets is computed for triangle of sizes h, I2, 13 
satisfying the conditions: 



li±ia<o.i ^i±ii<o.i, ^<o.i 



(15) 



rs ri 7-2 

Each triangle is ordered such that 7-3 < r2 < ri. We measure 
Fi for triangles with r-j between 5 arcsec and 80 arcmin, ranging 
in logarithmic bins of width log(r3) = 0.1. With this choice, we 
can assume that the distance between each pair of galaxies within 
the boxes that we are correlating has a maximum error of one bin. 
Finally we integrate Fi using equation (45) and equation (50) of 
Jarvis et al. (2004). As we already anticipated, the filter defined by 
Equation (13) for a given size & has infinite support so we ideally 
would need to measure Fi at all scales to perform the integral giv- 
ing the components of the the third-order moment of the aperture 
mass. However the integral is significant only for triangles up to 
r3 ~ 4i?, implying that our measurement of the third order moment 
is reliable up to an angular scale 1? ~ 20 arcmin. We tested our al- 
gorithm against the direct measurement of the third-order moment 
of the aperture mass on ACDM simulations in fields of 25 deg^ 



and we found indeed a good agreement between the two methods 
at these scales. 

Throughout the paper we use only the aperture mass statistics 
and we call the systematic produced by the intrinsic alignments on 
the second and third order moments of the aperture mass II and 
III, respectively. Similarly, we call GI the second order moment 
measured by using the filter defined by the Equation (13) produced 
by the shear-shape coupUng and we call the third order moments 
GGI and Gil. Finally, we call the second and third order moments 
of aperture mass produced by a weak lensing field GO and GOG. 

In order to know the importance of the GGI, GII and III when 
estimating the third order moment of the aperture mass GGG using 
the observed ellipticity of galaxies we have to compare those terms 
with the third order moment of the aperture mass produced by a 
weak lensing shear field. In the quasi-Unear perturbation theory, 
the perturbation of the density field S is considered to be small so 
that it can be developed in a series 5 = S^^^ + 5^^^ + where 
5^^^ is the linear evolving density field and 5'"' oc In 
this approximation the (Mfp)^ weak lensing signal is (Fry 1984, 
Schneider et al. 1998) : 



GGG = (Mfp)^ 



1 81Hg^3 
2n 4c6 " 



dw- 



poo / Q \ 

Jo ^ 

poo / Q \ 

Jo / 

X/(|S1 -|-S2|7?)F2(S1,S2), 



(16) 



with 



l-WH 

g{w) = / dw'ps{w' 

J w 



f^{w - w') 



(17) 



is the comoving angular diameter distance, P{s,w) is the 
3 dimensional power spectrum of matter fluctuations, Ps{w) is the 
comoving distance distribution of the sources, /(s??) is the Fourier 
transform of the filter as defined by Equation (14) and F2(si, S2) 
is the coupling between two different modes of density fluctua- 
tions characterised by the wave vectors si and S2. We compute 
^2(81,82) using the fitting formula suggested by Scoccimarro & 
Couchman (2001). 



3 N-BODY SIMULATIONS AND GALAXY MODELS 

In this analysis we use the same set of simulations used as in H06. 
We recall here only the main characteristics and refer the reader to 
H06 for more detailed information. The N-body simulation is a box 
of 300/r^Mpc reahzed using a ACDM cosmology. The matter 
density parameter is Qm = 0.3, the cosmological constant is Qa = 
0.7, the normalisation of the power spectrum of matter fluctuations 
is as = 0.8, the reduced Hubble constant is h = 0.7 and the 
baryon density parameter is fibft^ = 0.2. The simulation started 
at redshift z — 60 and evolved to 2: = 0. Twelve lines of sight 
were built by stacking boxes back along the 2-axis with random 
shifts between the boxes in the x and y direction in order to avoid 
artificial correlations, producing (almost) independent realisations. 

The over-densities are identified using a 'friends of friends' 
group finder, which allows one to identify the halos. The particle 
mass of 1.7 x lO'^" h~^MQ allows us to find bound halos with 
masses larger than few 10^^ h~^MQ. The halos are then populated 
with galaxies. The luminosity is assigned following the conditional 
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Figure 1. Left panel: (M|p ) signal generated by the intrinsic alignments (III) as a function of the angular size By comparing the E-modes (black diamond) 
with the expected shear signal (black solid line) one can determine the level of contamination produced by the existence of intrinsic alignments between 
galaxies. The shear signal has been computed using Equation (16) for the same redshift distribution used to compute the intrinsic alignment term III. The 
median redshift of this distribution is ~ 0.7. The error-bars are the dispersion of the average signal measured in twelve independent realizations. The 
B-mode signal, (M^ ), and the E/B components, (MfpMx ) and (MapM^ }, are shown dashed and for clarity we chose not to plot the error bars. The right 
panel shows the same results as the left panel for a lower redshift sample of galaxies < z < 0.4. For this redshift distribution characterised by — 0.3, 
the predicted shear signal GGG ~ 10"^ is indistinguishable from zero in this plot. 



luminosity function (CLF) of Cooray & Milosavljevic (2005) and 
the ellipticity is assigned using a toy model that has been shown 
effective at reproducing the observations of Mandelbaum et al. 
(2006) and Hirata et al. (2007). In this model, elliptical galaxies 
are given the same ellipticity as their parent halos. Spiral galaxies 
are modeled as a thick disk oriented almost perpendicular to the an- 
gular momentum vector of the halo with a mean random misalign- 
ment of ^ 20 deg (van den Bosch et al. 2002, Heymans et al. 2004). 



At each identified dark matter halo in the simulations we generate 
both a spiral and an elliptical galaxy so that we can investigate the 
results as a function of morphology. We define the ellipticity of the 
model galaxies using Equation (2). The resulting ellipticity distri- 
bution has a zero mean and dispersion of CTe = 0.33. 

In the results that follow we present three models: one com- 
posed exclusively of spiral galaxies, one composed exclusively of 
elliptical galaxies and one containing both the morphologies de- 
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noted 'mixed'. The mixed model has been buiU following Cooray 
& Milosavljevic (2005) and it contains ~ 30% elliptical galaxies. 
The mixed model is the most realistic model out of the three tested; 
the GI and II components predicted using the mixed model agree 
with the ones measured using the SDSS data (Mandelbaum et al. 
2006, Heymans et al. 2006, Hirata et al. 2007). The number density 
of the final catalogue, which contains galaxies between < z < 1 
and is complete up to r = 25.5, is about 5 per arcmin^, which is 
considerably lower than what is expected for such a survey. This 
is due on one hand to the fact that the low resolution of the sim- 
ulations implies the loss of low-mass halos (only halos with mass 
larger than few lO^'^ M© are identified) and that each halo is pop- 
ulated with only one galaxy. The lack of the low-mass halos and 
satellite galaxies is the main limitation of our model and this point 
is discussed further in the conclusions. 

Finally, these same N-body simulation which have been popu- 
lated with galaxies are ray-traced to produce twelve projected mass 
distributions n for a source plane at Zs = 0.45 and for Zg — 1.05. 
We use these simulated k maps in section 5 to study the shear-shape 
effect. Each projected mass map covers an area of 5 x 5 deg^ in 
2048 X 2048 pixels. 



4 INTRINSIC ELLIPTICITY ALIGNMENT 

In this section we report the third order moment components of the 
aperture mass given by the intrinsic alignment of the sources. We 
measure the third order moment of the aperture mass for a survey 
with maximum depth z ~ 1 and compare the measured intrinsic 
alignment to the expected weak lensing GGG signal from Equation 
(16). The redshift distribution of the sources Psiw) is computed 
directly from the catalogues used for the III measurement so that 
both III and GGG terms have been computed for the same survey. 

Figure 1 shows the three-point intrinsic alignment for our 
three toy models: elliptical (upper panel), spiral (middle panel) and 
mixed (lower panel) for two survey types; one with a median red- 
shift of «m ~ 0.7 (left panels) and the other with Zm 0.3 (right 
panels). The error-bars on the measurements have been computed 
as the dispersion in the average signal measured in the twelve real- 
isations, thus they include statistical and cosmic variance. For the 
moderately deep survey (zm ~ 0.7) the III term is consistent with 
zero for all the three models. This result is in agreement with H06 
who find that the two-point intrinsic ellipticity term II, is consistent 
with zero for the same survey characteristics. As expected, how- 
ever, we find that when the redshift depth of the survey decreases 
the III term becomes significant due to the fact that for lower red- 
shift galaxies, a given angular distance corresponds to a physically 
closer triplet. Comparing now the results for the three different toy 
models for this shallow survey we find very different results. Fig- 
ure 1 shows that the III term is significantly positive for i9 > 2 
arcmin for elliptical galaxies (upper right panel) meanwhile it is 
slightly negative for the spiral model (middle right panel) and sig- 
nificantly negative for the mixed model (lower right panel). The 
mixed model contains roughly 70% of spiral galaxies and 30% of 
elliptical galaxies, such that the most frequent triangle in the three- 
point measurement contains two spiral galaxies and one elliptical 
galaxy. As this type of correlation does not exist when one consid- 
ers only elliptical or only spiral morphologies, the mixed result is 
not a weighted average of the result containing only elliptical or 
spiral galaxies. We verified that the third order moment of the aper- 
ture mass is indeed negative when we correlate triplets containing 
two spirals and one elliptical. Figure I shows that for shallow sur- 
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Figure 2. Left side: ratio between the intrinsic alignment (II) and the aper- 
ture mass variance (GG) for three surveys. Right side: ratio between III the 
expected weak lensing signal (GGG) for the same surveys as the left panel. 
The galaxy model is a mixed model. The expected second and third order 
aperture mass moment are computed for each survey by using the same 
source redshift distribution used to measure the III and II terms. These dis- 
tributions are characterized by: Zm ~ 0.31 for source distribution between 
< z < 0.4, ~ 0.4 for source distribution < z < 0.5 and source 
distribution Zm ~ 0.6 for < 2 < 0.8. 

veys (Zm ^0.5) the intrinsic alignment dominates the signal and 
suggests that also for deeper surveys the intrinsic alignment could 
affect the measured the three-point weak lensing signal. 

In order to quantify the effect of the intrinsic alignment on 
the weak lensing shear statistics we compare the amplitude of the 
two- and three-point intrinsic alignment signal II and III with the 
expected weak lensing signal, respectively GG and GGG, for dif- 
ferent shallow surveys. Figure 2 shows the amplitude of the II/GG 
(left panel) and III/GGG (right panel) ratios for three different sur- 
veys depth using the mixed model as an example. We find that for 
a given survey depth, the ratio between the III term and expected 
weak lensing third order moment GGG, is generally higher than the 
ratio between the II term and the variance of the weak lensing aper- 
ture mass GG. One can think to reduce significatively the intrinsic 
ellipticity contamination to the two- and three-point shear statistics 
by removing very low redshift sources (z ;S 0.2 ) which contribute 
to the intrinsic alignment signal but not to the weak lensing signal. 
We checked that the ratio II/GG and III/GGG drops on small scales 
1? ^ 10 arcmin if galaxies with z < 0.2 are rejected, but it is almost 
unchanged at larger scales. 

Figure 1 and Figure 2 show that without a technique to cor- 
rect for the intrinsic alignment it will be impossible have a precise 
measurement of the three-point shear statistics in shallow surveys 
(zm 0.5). For moderately deep surveys like the CFHTLS Wide, 
the ratio III/GGG is consistent with zero as shown by the left panels 
of Figure 1. However, for larger surveys, where the statistical and 
cosmic variance are small the intrinsic alignment could still play a 
role in the precision with which one can constrain the cosmologi- 
cal parameters. This can be seen in Figure 2 where for a relatively 
deep survey (z < 0.8) the intrinsic alignment contribution is likely 
to be non-zero, specially for the third order moment of the aperture 
mass. 

This result strongly supports the need for reliable photometric 
redshift for weak lensing studies. Moreover, Figure 2 demonstrates 
that if the redshift is known, one can select sources in order to en- 
hance the II and III signal relative to the contribution of the GG 
and GGG terms to the measured two and three-point shear statis- 
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tics. Tiiis is an important result since one needs a good intrinsic 
alignment model to be able to remove this contribution from the 
two-point (three-point) weak lensing shear statistics (King 2005). 
The fact that the in/GGG signal is more significant than the II/GG 
signal implies that it may be easier to study intrinsic alignments us- 
ing three-point statistics. In this respect Figure 2 shows that large 
multi-band surveys like SDSS and KIDS or PanSTARRS, represent 
excellent surveys to investigate and model the intrinsic aUgimient 
of elliptical and spiral galaxies. 



5 SHEAR ELLIPTICITY ALIGNMENT 

In this section we study the behavior of the three-point shear-shape 
coupling which, as it has been detailed in equations (6) and (7), can 
be divided in two terms, namely GGl and Gil. One would expect 
these two terms to have a different behavior as a function of the 
redshift distribution of the lenses and of the sources. The GII term 
depends on the intrinsic ellipticity correlation between two pairs of 
galaxies thus it should be significant only for triplets which contain 
foreground galaxies closer than the scales on which tidal forces act. 
The GGI term correlates the shape of lensed backgrovmd galaxies 
with the intrinsic shape of a foregrovmd galaxy so it is expected to 
be significant also when correlating well-separated redshift slices, 
similarly to what has been found for the GI term (Bridle & King 
2007). Using both ray tracing simulations with redshift Zs = 0.45 
and 2s = 1-05 and dividing the foreground galaxies into redshift 
bins, we study the evolution of the GII and GGI terms as a func- 
tion of the source and lens redshift distribution. Figure 3 shows 
the amplitude of the GGI (left panels) and Gil (right panels) terms 
for elliptical (upper panels), spiral (middle panels) and the mixed 
model (lower panels) for different redshift bins. We note that the 
amplitude of the GGI and GII terms for a given survey depends 
on the morphology of the galaxies. The elliptical galaxies show a 
significant negative GGI signal whereas the spiral model shows a 
GGI signal which is consistent with zero. Finally, the mixed model 
shows a slightly negative GGI signal. 

The GII component for the spiral galaxies show an angular 
dependence similar to the one of the elliptical galaxies and inter- 
estingly the signal at intermediates scales (2 <, i9 <, 10 arcmin) is 
stronger than for the elliptical galaxy model. Similarly to what we 
previously found for the III term, the mixed model GII term is dif- 
ferent from the one obtained using only elliptical or spiral galaxies. 
This is likely to be a consequence of the fact that for the two mor- 
phologies the ellipticity depends differently on the properties of the 
parent halos, which also could explain the difference in results for 
the III term. To understand this phenomena would require build- 
ing a model for the three-point correlation function between weak 
lensing shear and tidal field, in a similar analysis as that of the two- 
point correlation function presented in Hirata & Seljak (2004). This 
however is beyond the scope of this paper. 

In order to allow a more quantitative comparison between the 
weak lensing three-point shear signal and the systematics Gil and 
GGI, we report in Figure 4 the ratio between the terms Gil and 
GGG (left panel) and GGI and GGG (right panel) for the mixed 
model. To aid future comparisons we provide fits to our results. We 
make the assumption that for a given triplet the GGI and GII terms 
can be factorized in two functions: one depending on the comoving 
distances of lenses and sources modeled as in King (2005) and the 
other depending on the angular scale. We rewrite the GGI term as: 



with: 

EgGi(Wsi, Wsa^WL) = / 

Jo 



Pi(wi), 



(19) 



/fc(WB2) 

where wl is the maximal comoving distance of the lenses with the 
condition wl < min(wsi, Wsj), fk{w) is the angular diameter 
distance, pi(wi) is the radial lens distribution and -F(i9) is a generic 
function of the angular scale. 

Similarly we factorize GII with the assumption that the intrin- 
sic aUgnment acts only for pairs of galaxies with the same redshift. 
Thus we write: 



GII(ws, WL, = Egu{w„wi^)F{'&), 



(20) 



with: 



dwi — — ^ pi(iui) , (21) 



with the condition Wh < Ws- For a broad source distribution one 
should integrate Equation (19) and Equation (21) over the source 
distribution to obtain the average effect. However in our case the 
ray-tracing, source planes follow a Dirac distribution centered at 
Zs = 0.45 or Zs = 1.05. The lens distribution pi(wi) is measured 
directly from the galaxy catalogues. We find that these scaling fac- 
tors provide a good description of the redshift dependence both for 
the GII and GGI term. With this redshift dependence model we try 
to model the behavior of the GII and GGI terms as a function of the 
angular scale. We use a two parameter function F(i?): 



F{d) = A exp(-t?/i?o) 



(22) 



GGI(Wbi, Wb2, = EgGi(Wsi, Ws2,ML)i^(^?), 



(18) 



where A and t?o are free parameters. 

In order to avoid bias from the limiting resolution of the simu- 
lations we chose to perform the fit using only angular scales i? > 1 
arcmin. In Figure 3 we compare the measurements of the GGI and 
Gn terms with the best-fit model (solid lines) obtained by using the 
redshift rescaling defined by the equations (19) and (21) and an an- 
gular scale dependence defined by Equation (22) for the elliptical, 
spiral and mixed galaxies. Figure 5 shows the best fit parameters A 
and t?o for the GGI (left panel) and the GII (right panel) terms for 
several redshift slices, for elliptical, spiral and mixed model. These 
values are also summarized in the tabic 1 which includes also the 
reduced of the fit for both the GGI and GII terms. The small 
values of the show that the model we suggest is a good fit to 
the measured GGI component. However, for the elliptical model 
the dispersion of the best-fit A parameter between the four red- 
shift slices (see left panel of Figure 5) is larger than the error bars. 
This suggests that the redshift-dependent rescaling described by the 
Equation (19) could require some modification. The GII compo- 
nent is generally more noisy than the GGI component, and for this 
reason it is hard to establish whether the model defined by equa- 
tions (21) and (22) is a good one. However, for the mixed model, 
which is our most realistic model, the GII term can be fairly well 
described by our fitting function. Increasing the size of simulations 
in future analyses will allow us to improve upon these models. 

Because of the different dependence on the morphology of 
galaxies and on the redshift distribution of the sources and lenses 
for the III, GGI and GII terms, one may be interested in knowing 
the total effect of the intrinsic alignments and of the shear-coupling 
on the three-point shear for realistic surveys. Figure 6, shows the 
ratio between the sum of GGI, GII and III terms and the expected 
weak lensing signal for several redshift distributions for both the 
elliptical (left panel) and mixed (right panel) model. 
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Figure 3. GGI (left panels) and GII (right panels) components for different redshift slices and three galaxy models: elliptical (upper panels), spiral (middle 
panels) and mixed (lower panel). We show three sets of measurements corresponding to: a lens distiibution < < 1 and a source distribution Zb = 1.05 
(black diamonds), < zi < 0.4 and Zg = 1.05 (red triangles) and < < 0.4 and Zg = 0.45 (blue squares). For the GGI component we also include 
a third model which contains lenses < < 0.4 and sources in two different planes: one at Zg = 0.45 and the other at Zb = 1.05. Error bars represent 
the dispersion in the average value between the twelve simulations. For each measurement of the GGI and GII components we show the best-fit model (solid 
lines) obtained by rescaling the redshift distribution dependence given by equations (19) and (21) with the angular dependence described by Equation (22). 



For shallow surveys (zi £^ 0.5), the ratio is largely dominated 
by the III term. This is true both for elliptical and mixed model 
with the difference that for the elliptical galaxies the intrinsic align- 
ment enhances the weak lensing shear signal whereas for the mixed 
model it suppresses the weak lensing signal. For a deep survey 
Zm — 0.7 the ratio is dominated by the GGI term, whose amplitude 
is 10% of the GGG signal for elliptical galaxies and few percent 
for mixed model. 

In order to compare our results to observations we present in 



in Figure 7 the {M'ip) results by Jarvis et al. (2004) and the (A/fp) 
that we would expect for this same survey from weak gravitational 
lensing and the III, GGI and GII terms. For this comparison we 
have used all of the three galaxy models. We compute the GGG 
signal using a redshift plane at Zs = 0.66 as done by Jarvis et al. 
(2004). We add the III model using galaxies with z < 0.8 so that 
the redshift distribution is characterized by the a median redshift 
similar to the one of the CTIO galaxy catalogue. We use Equation 
(19) and Equation (21) to rescale the GII and GGI to a survey with 
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Table 1. Summary table of the best fit parameters of the GGI (first four columns) and GII (last three columns) terms for elliptical (upper lines), spiral (middle 
lines) and mixed (lower lines) galaxies using Equation (19) and Equation (21). The values the parameters i9o are given in arcmin, the values of the parameter 
A are given in units of I0~^h Mpc~^. For each model we report the value of the reduced x^. 
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Figure 4. Left side: ratio between the GII term and the expected signal 
(GGG) for three shallow surveys. Right side: ratio between the GGI term 
and the GGG signal for the same surveys as the left panel. The galaxy model 
is a mixed model. We used the same surveys as Figure 3. The GGG signal 
has been computed by taking Zs = 1.05 (black diamonds and red triangles) 
and 2s = 0.45 (blue squares). 

Zs ~ 0.66 and zi < 0.66. For this survey we find that the dominant 
contribution is given by the GGI which slightly decreases the esti- 
mated the weak lensing signal. We find that all three galaxy models 
are consistent with the Jarvis et al. (2004) results. 



6 CONCLUSION 

Using a set of realistic N-body ACDM simulations we have ex- 
plored the effect of intrinsic galaxy alignments and the coupling 
between weak lensing and the foreground ellipticity field on the 
third order moment of the aperture mass. We find that the intrin- 
sic alignment dominates the three-point shear signal for shallow 
surveys. For deeper surveys the intrinsic alignment is less signif- 
icant, as a result of projection effects. Nevertheless, if not taken 
into account properly, it will still limit the accuracy of tomography 
measurement and affect the predictions of cosmological values for 
the next generation of very large surveys, where the signal to noise 
ratio is high. 

We found that for a given survey depth intrinsic alignments 



affect the three-point weak lensing statistics more strongly than 
the two-point shear statistics. In other words, in order to achieve 
the same level of accuracy in the three-point shear statistics mea- 
surement as the two-point weak lensing shear statistics one needs 
deeper surveys. 

Overall, this result shows once more the importance of the 
knowledge of redshift for each individual galaxy for high-precision 
cosmology. Knowing the redshift of each source allows one to re- 
move the bias from intrinsic galaxy alignments by discarding phys- 
ically close triplets (pairs) when computing the three-point (two- 
point) shear statistics. Moreover the knowledge of redshift allows 
one to model the intrinsic alignment and remove its effect on the 
two-point weak lensing shear statistics (Joachimi & Schneider, in 
prep.). 

In this perspective, the measurement of the three-point shear 
statistics from a shallow survey offers a particularly effective way 
to test intrinsic alignment models. As we showed in this work, 
the three-point shear statistics in low redshift bins are dominated 
by the intrinsic alignment term, permitting accurate measurements 
on intrinsic alignment models. As future dark energy surveys will 
rely heavily on the good modeling of these effects so they can be 
marginalised out, this result provides an important new route to 
constrain and model this physical systematic effect. It is indeed 
possible to choose low-redshift bins in order to enhance the intrin- 
sic alignment signal so that the shear signal becomes negligible; for 
example, by selecting galaxies at z < 0.4. In this case the III/GGG 
ratio is around fifty, allowing one to study the intrinsic alignment, 
essentially without contamination from the weak lensing shear. 

Using projected mass maps we have studied the shear-shape 
coupling effect which is also likely to bias the measurement of the 
third order moment of weak lensing shear. We showed that this 
systematic can be described by two terms, GGI and GII, which af- 
fects the third order moment measurement in a non-trivial way and 
is dependent on the distance between the lenses and sources and 
on the morphology of galaxies. Even for a moderately deep survey 
like the CFHTLS Wide the amplitude of the third order moment 
of the shear estimation could be underestimated by ~ 5 — 10%. 
For shallower surveys such as KIDS or PanSTARRS-1 the bias is 
expected to be higher. Our results show that it will not be possible 
to carry out precise three-point cosmic shear measurements with 
these surveys without modeling the coupling between weak lens- 
ing shear and intrinsic alignment. Also for the next generation of 
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Figure 5. Left panel: best-fit values for the parameters i9o and A of Equation (22) used to fit the angular dependence of GGI for the four redshift distributions; 
we chose to indicate each source/redshift distribution using the same symbols as Figure 3. Right panel: shows the same parameters as left panel but now for 
the GII component. 



o 
o 



10 

0.1 ^ 
-0.1 

-10 



elliptical 



^ ^ ^ ^ 



10. 



'i)'(arcmin) 



o 
o 
o 



o 
o 



10 ^ 
0.1 

-0.1 
-10 t 



1. 



mixed 



Oz,= 1.0, 0<Z|<1.0 
Az,= 1.0, 0<z,<0,4 
□ z =0.5, 0<z,<0.4 




10. 



i?(arcmin) 



Figure 6. Ratio between the sum of the terms Oil, GGI and III and the weak lensing third order moment signal GGG, for the elliptical (left panel) and mixed 
model (right panel). Three cases are shown: = 1, < zi < 1 (black diamonds), Zs = 1, < < 0.4 (red triangles) and Zb = 0.5, < Zi < 0.4 (blue 
squares). For the first and second case the GGG component is computed assuming Zs = 1.05. For the third case Zs = 0.45. 



deep large surveys, such as SNAP, DUNE or LSST the precision 
one can achieve on the cosmological constraints relies in the ability 
to model and marginalize out the shear-shape and intrinsic correla- 
tions. 

It is possible to use simple models to describe both the angular 
and redshift dependence of the GGI and Gil terms. These models 
can be used to correct the effect of the coupling between shear and 
intrinsic alignment on the third order moment of the aperture mass. 

Unfortunately, the sample available for this work is not big 
enough to give a definitive answer; more precisely, the parametric 
models we used to fit the GGI and Gil components are marginally 
constrained, due to the large statistical and sampling variance af- 
fecting our sample. A more detailed study of the redshift depen- 
dence of the shear-ellipticity correlation will be required. Such a 
study should use a larger sample of simulations in order to bet- 
ter constrain our models. Furthermore, we still need to develop a 
method to include satellite galaxies and account for environmen- 
tal effects when determining the ellipticity of galaxies in a single 
parent halo. In addition, an improved mixed galaxy population as a 
function of redshift evolution would also play an important role in 
establishing the correct average effects, as we have shown in this 
paper the net effect depends on the morphology of the galaxy pop- 
ulation. We also think that a comparison between results from sim- 
ulations with results from real data is the best way to validate our 



results. Concerning this point, the agreement between simulations 
and real data shown in the study of the 11 (H06) and Gl (Hirata et 
al. 2007) components is a good indication of the fact that the sim- 
ulation used in this work, even if it represents a simplified model, 
shows the correct dependence of the intrinsic alignment effect on 
morphology of the galaxies, luminosity and redshift. 

Like other previous works on intrinsic alignment and shear- 
shape coupling, this paper demonstrates the great importance of re- 
liable redshifts for future and current galaxies surveys, which could 
be used, in parallel with simulated catalogues, to study and correct 
for intrinsic alignment and shear-shape coupling on the two- and 
three-point weak lensing statistics. 
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